Interior design of a two-dimensional semiclassic black hole 



Dana Levanony and Amos Ori 
Department of Physics, Technion, Haifa 32000, Israel 
(Dated: December 28, 2009) 

We look into the inner structure of a two-dimensional dilatonic evaporating black hole. We 
establish and employ the homogenous approximation for the black-hole interior. The field equations 
admit two types of singularities, and their local asymptotic structure is investigated. One of these 
singularities is found to develop, as a spacelike singularity, inside the black hole. We then study the 
internal structure of the evaporating black hole from the horizon to the singularity. 
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INTRODUCTION 

The discovery of Hawking radiation and the black- 
hole (BH) evaporation phenomenon [l[ raised several out- 
standing open questions and intriguing puzzles. One such 
problem which attracted a lot of attention is the infor- 
mation puzzle: Simple thought experiments suggest that 
in the process of black-hole formation and evaporation, a 
pure initial quantum state will evolve into a mixed state, 
and consequently part of the information encoded in the 
initial state will get lost. 

Apparently this problem of information loss is inti- 
mately related to another conceptual problem in black- 
hole physics: The formation of a spacetime singularity 
inside the BH. Pictorially speaking, a portion of the ini- 
tial information propagates to the singularity and dis- 
appears there. It is widely anticipated, however, that 
the formation of spacetime singularities is a mere artifact 
of the classical (and semiclassical) theory, but Quantum 
Gravity will eventually resolve the black-hole singulari- 
ties, and presumably solve the information puzzle as well. 

In 1992, Callan, Giddings, Harvey, and Stromingcr 
(CGHS) [2j proposed a simplified framework for inves- 
tigating black-hole formation and evaporation. They in- 
troduced a two-dimensional toy-model with gravity and 
matter fields coupled to a dilaton scalar field. At the 
classical level this model admits a one-parameter fam- 
ily of static black-hole solutions (parameterized by their 
mass). When semiclassical corrections are added, the 
two-dimensional black hole emits thermal radiation and 
consequently evaporates. Based on the trace anomaly, 
CGHS provided simple explicit expressions for the semi- 
classical contribution to the energy-momentum tensor. 
Consequence the semiclassical dynamics may be formu- 
lated as a closed (constrained) system of second-order 
partial differential equations. It thus provides a sim- 
ple framework for exploring various aspects of black-hole 
evaporation, and particularly the information puzzle Q • 

Originally it was hoped Q that the CGHS evaporat- 
ing BH will be free of any singularities. However, Russo, 
Susskind and Thorlacius [H soon found that a singular- 
ity inevitably develops inside the CGHS BH, at a certain 
value of the dilaton filed. Thus, the CGHS formalism 



does not resolve the information puzzle at the semiclas- 
sical level. One may still hope, however, that when the 
model is fully quantized the singularity (and the informa- 
tion puzzle) will be resolved. This approach was pursued 
by Ashtekar, Taveras, and Varadarajan who formu- 
lated the quantum-field analog of the two-dimensional 
CGHS model. In this set-up, the original CGHS variables 
are replaced by quantum operators. Spacetime evolution 
is then described by a system of operator partial differ- 
ential equations. The exact solution to these operator 
equations is not known. Nevertheless, some approximate 
solutions were constructed [f|, illuminating certain as- 
pects of the problem of black-hole evaporation, and sug- 
gesting that the singularity will indeed be resolved in the 
quantized theory. 

The above discussion signifies the inevitable spacetime 
singularity as a key feature of the CGHS semiclassical BH 
spacetime. However, the detailed structure and prop- 
erties of this singularity have not been investigated so 
far to the best of our knowledge. The main purpose 
of this manuscript is to present such a detailed analysis 
of the structure of that singularity. The motivation for 
this investigation is obvious from the discussion above: 
First, since this singularity plays such a crucial role in 
the information-loss puzzle, it will be useful to better 
understand its properties. For example, one would like 
to know how strong this singularity is, and what are the 
prospects for extending semiclassical spacetime beyond 
it. Second, understanding the asymptotic behavior of 
the various fields at the singularity may provide a use- 
ful starting point for exploring how quantum treatments 
(like that of Ref. [B| for example) may resolve the singu- 
larity. 

Although the main objective of this paper is the 
asymptotic behavior near the singularity, we also ana- 
lyze here the internal structure of the evaporating semi- 
classical BH in the entire range from the horizon to the 
singularity. We do this by constructing approximate so- 
lutions in various domains of the BH interior, and then 
matching these solutions at their respective overlap re- 
gions. Understanding the entire BH internal structure 
is interesting on its own right, but is also important for 
full determination of the singularity structure: When the 
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latter is derived by a purely local analysis, one obtains 
a family of local asymptotic solutions which depend on 
certain free parameters (or free functions, depending on 
the context). The value of these parameters in an actual 
BH solution needs to be determined by matching to ini- 
tial conditions (e.g. at the horizon). The full internal 
solution, from the horizon to the singularity, is required 
for obtaining the right values of these parameters. 

As long as the evaporating black hole is macroscopic, 
the semiclassical effects are very weak in a local dynam- 
ical sense. This allows one to employ the adiabatic ap- 
proximation. Namely, any local region of spacetime may 
be well approximated by a certain classical CGHS BH so- 
lution (with a certain mass parameter). The classical BH 
solutions are static outside the horizon but homogeneous 
inside it, implying that the interior of a macroscopic 
evaporating BH is (locally) approximately homogeneous. 
This homogeneous approximation is a key element in our 
investigation. It greatly simplifies the analysis, as the 
field equations now reduce to ordinary differential equa- 
tions. 

In Sec. II we present the CGHS action and field equa- 
tions, and re-formulate them in new variables. Sec. Ill is 
devoted to the homogenous approximation. We explore 
its main properties and its domain of validity. In Sec. 
IV we analyze the singularities of the homogenous field 
equations, and obtain the local asymptotic behavior of 
the solution near these singularities. In Sec. V we con- 
struct an approximate solution for the fields inside the 
evaporating BH, in the entire range from the horizon to 
the spacelike singularity. 

In Sec. VI we present simple effective Lagrangian and 
Hamiltonian which yields the local asymptotic dynam- 
ics near the spacelike singularity. The construction of 
such effective Lagrangian and Hamiltonian is motivated 
by the idea, that perhaps it will be possible to employ 
simple quantum-mechanical considerations to resolve the 
singularity (and thereby to explore the possible exten- 
sion of spacetime beyond it). Although such a quantum- 
mechanical treatment is beyond the scope of the present 
paper, we nevertheless take here the preparatory step 
(carried out entirely within the semiclassical framework) 
of constructing the effective near-singularity Lagrangian 
and Hamiltonian. Finally, in the last section we summa- 
rize our main results and discuss their significance. 



THE MODEL AND FIELD EQUATIONS 

We begin by presenting the two-dimensional semiclas- 
sical CGHS model Q. Beside the dilaton </>, this model 
contains a large number N > > 1 of identical scalar mat- 
ter fields fi , and a cosmological constant A 2 . We express 
the metric in double- null coordinates u,v (the "confor- 
mal gauge"), namely ds 2 = —e 2p dudv. The action then 
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The last term in the action expresses the semiclassical 
effects. The Einstein equations take the form: 

T uv = e~ 2 * (2(j), uv -4<f>, u <f>„ ~\ 2 e 2 P) - £p mv = (2) 

and 

1 N 

Tun = e 20 (4p, u cj), u -2<j>, uu ) + - 2J fi,u fi,u (3) 

i=l 

N 

~^}P™ P™ ~Pmu -z u (u)\ = 0, 

1 N 

T vv = e 20 (4p, v <f>, v -2(f> rvv ) + - 2J fi, v fi,v (4) 

i=l 

N 

-J2 \- P}V Piv ~ p,vv = °> 

where z u (it) and z v (v) encode the information about the 
initial quantum state, which in turn determines the semi- 
classical fluxes. The dilaton and matter equations are: 

- 4(f>, uv +4<p, u (f>, v +2p mv +X 2 e 2 P = 0, (5) 
dudvfi = 0. 

Throughout this paper we set A = 1. This choice is 
equivalent to the change of variable p' = p + ln(A), which 
does not affect the field equations otherwise. We also 
set fi — 0, as we are dealing here with the evaporation 
rather than formation of the BH. 

Following Ref. Q, we define new variables: R = e~ 20 , 
S = 2{fi — <f>). A straightforward substitution of the new 
variables in Eqs. (|2I5|) yields: 
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where p = h(S— In R) is to be substituted, and K = N/12 
expresses the magnitude of the quantum effects. The 
constraint equations become: 



R. 
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where hereafter w stands for either u ov v. The semi- 
classical energy fluxes T ww along both null directions are 
given by: 
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It is useful to re-express the system of evolution equa- 
tions (J6j) in its standard form, in which R ruv and S ruv are 
explicitly given in terms of lower-order derivatives: 

-n-jmi — 2(R-K) 2R(R-K) ' W 



9 - e S 



K 

2R(R-K) 



R,u R,v 2R 2 {R-K) ■ 



3 



This form makes it obvious that the evolution equations 
become singular when R = K . This singularity, which in 
the original variables takes place at <j> = — |ln(^|), was 
already noticed preivously Below we shall explore in 
some detail the homogeneous variant of this singularity. 

The semi-classical equations reduce to the classical 
ones by setting K = 0: 

R,uv = ~e s , (10) 
<? — n 

and the constraint equations: 

This set of equations admits a one-parameter family 
of solutions (up to gauge transformations), which is 
the dilatonic two-dimensional analog of the standard 
Schwarzschild solution. Each member of this family 
describes a black-hole space-time, which (like its stan- 
dard four-dimensional counterpart) is static outside the 
BH and homogeneous inside it. Throughout this paper 
we shall simply refer to this class of solutions as the 
" Schwarzschild solution" (despite a slight abuse of stan- 
dard terminology). We shall focus on the internal part 
of the BH. In Eddington-likc double-null coordinates the 
internal solution takes the form 

R = -e v+u + M, (11) 
S = v + u, 

where M is the (dilatonic, 2-dimensional) Schwarzschild 
mass parameter. 

HOMOGENOUS SET-UP 
Justification and domain of validity 

We shall now restrict our attention to the homogeneous 
solutions of the semiclassical field equations inside the 
evaporating BH, namely, solutions which only depend on 
t, where we define 

t = v + u, x = v — u. 

We first need to discuss the justification for this homo- 
geneous approximation and its domain of validity. The 
classical interior solution (fTTj) is obviously homogeneous 
as it depends on v + u solely. In an evaporating BH, 
however, the semiclassical effects spoil the exact homo- 
geneity, as manifested by the drift in the BH mass. Yet, 
as long as the BH is macroscopic (M >> K), this drift 
is very slow. Indeed it undoubtedly has a dramatic effect 
over the long evaporation time-scale — the BH eventually 
disappears (or at least becomes microscopic) after all. 
Yet, from a local point of view this drift is of negligibly 



small rate. Thus, over a typical dynamical time/length 
scale (say, Av,Au of order ~ 10), the change in M is 
SM << M and may be neglected, which naturally leads 
to the homogeneous approximation. 

The local, relative magnitude of the semiclassical terms 
is of order K/R, as can be seen in the evolution equations 
(|9|). As long as the BH is macroscopic (M >> K), and 
as long as we are dealing with the portion R >> K oi the 
BH space-time, we may apply the adiabatic point of view 
to the evaporation process: Namely, in each region of 
space-time the solution is approximately Schwarzschild, 
with a local effective mass parameter M which slowly 
drifts in v and/or u. We shall refer to this slowly varying 
function M(v,u) as the effective mass of the evaporat- 
ing BH. (In fact, as long as the interior of the BH is 
considered, the effective mass depends only on v. (llj ) 
The adiabatic approximation, wherever applicable, au- 
tomatically implies approximate homogeneity inside the 
BH — simply because the Schwarzschild interior solution 
is homogeneous. 

It should be emphasized, however, that the adiabatic 
approximation docs not apply in the small- R region near 
the singularity, where R becomes comparable to K. In 
this region there are strong local semiclassical effects, and 
the local solution is very different from Schwarzschild 
(as demonstrated in the next section). We argue, how- 
ever, that the homogeneous approximation is still valid 
in this small- R domain. Generally speaking, this follows 
from causality: Since the small- R region is located at the 
causal future of the moderate- R interior region, it simply 
inherits the approximate homogeneity of the latter. 

To be more specific, let us pick a point inside the evap- 
orating BH, in a region where the effective mass is still 
macroscopic, M{v) = Mq >> K. This point is picked 
at a certain R value, R = i?o, which is >> K — say 
i?o = Mo/2 (recall that the horizon is located at R « M). 
Let us denote the coordinates of this point by (vq,Uo), 
and let to = vq + uq.xq = vq — uq. Let £ denote an 
initial surface which is the line t = to restricted to the 
range x$ < x < xq + Ax, with Aa; taken to be ~ 10 
(say). This construction is illustrated in FigQ] We shall 
now consider two cases, namely two slightly different so- 
lutions: In case (i), we specify on E initial data which 
exactly correspond to the classical Schwarzschild solu- 
tion with mass M(vq) = Mq. Since these data are inde- 
pendent of x by construction, the evolving solution will 
be precisely homogeneous, throughout the domain of de- 
pendence D + (E). Provided that Ax is taken to be large 
enough, D + (T,) will reach the small- R space- like singu- 
larity (more precisely, the R = K singularity which we 
explore in the next section). To this end Aa; only need 
to exceed the span of t from to to the singularity (of 
the homogeneous solution) . This span is of order 1 (for 
instance, for i?o = -Mo/2 this t-span is Zn(l/2) in the 
classical solution, and approximately the same number 
in the semiclassical solution). Hence taking Ax ~ 10 
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(say) guaranties that Z?+(E) will include a portion of the 
small- R singularity (as illustrated in FigfT]). 




FIG. 1: Spacetime diagram of a CGHS black hole which forms 
by gravitational collapse and subsequently evaporates. The 
thick null line at the bottom left represents the collapsing 
massive shell. The event horizon is marked by a dotted line. 
The spacelike hypersurface E is centered at v — vq (where the 
effective mass is M(v) = Mo), and is taken to be at a certain 
R value satisfying K « R < Mo. The gray area represents 
D+(E). The point E marks the end-point of the evaporation 
process. 

In case (ii) , we specify on E the initial data which cor- 
respond to the actual evaporating BH solution. Since 
both Mq and Ro arc >> the adiabatic approxima- 
tion applies, and the solution in the neighborhood of 
(i>o,uo) is approximately Schwarzschild. Therefore, we 
may regard the actual initial data (ii) as the homoge- 
neous, Schwarzschild data of case (i) plus a small pertur- 
bation. This initial perturbation will evolve into a small 
perturbation throughout the interior of Z)+(E). The rel- 
ative magnitude of this perturbation will be controlled 
by the small parameter K/Mo, and for sufficiently small 
value of this parameter we may neglect the perturbation 
at the leading order. 

The smallness of the perturbations in the interior of 
D + (E) is then guaranteed by standard stability theorems 
for nonlinear hyperbolic systems. One may be concerned, 
however, about the effect of the perturbation in the im- 
mediate neighborhood of the small- R singularity, and in 
particular on the very structure of the latter. A close 
examination, which is beyond the scope of the present 
paper, reveals that this singularity is locally stable to in- 
homogeneities. (That is, an inhomogeneous variant of 



the R = K singularity exists [7| , and this inhomogeneous 
singularity is generic, in the sense that it depends on four 
arbitrary functions of x [13].) 

We conclude that in the macroscopic-mass domain 
(Mo >> K) the homogeneous solution provides a good 
approximation not only in the domain R » K, but also 
up to the small- R singularity. This was also verified nu- 
merically by directly integrating the field equations ([6]) 
with initial data corresponding to a BH which forms by a 
collapsing null shell and subsequently evaporates Q. In 
particular it was numerically verified that the structure 
of the evolving spacelike singularity is well described by 
the homogeneous R = K singularity (described in the 
next section). 

In the above construction we considered for concrete- 
ness an initial hypersurface E located at R = Mo/2. We 
could have started instead at any other R = Rq initial 
value which satisfied K « Rq < Mq. The evolving ho- 
mogeneous solution depends very weakly on Rq , because 
of the local smallness of the semiclassical corrections at 
R >> K. For concreteness, and in order to avoid the 
arbitrariness associated with the extra parameter Rq, in 
the analysis below we shall consider the homogeneous 
solution obtained at the limit Rq — ► Mo. The solution 
obtained in this way may be thought of as the (horizon- 
regular) semiclassical counterpart of the homogeneous in- 
terior Schwarzschild solution. [l3[ ft carries (besides K) 
a single parameter Ma, representing the BH's remaining 
mass at the epoch of interest. In the rest of the paper we 
shall simply refer to this parameter as M for brevity. 

Homogenous field equations 

In the homogeneous framework the evolution equations 
([9]) become 

^ L 2(R-K) 11 2R(R-K)' 

q — r S K I p2 K ho\ 

° — c 2R,(R~K) x 1X 2R 2 (R-K) ' ^ ±J J 

where an over-dot denotes differentiation with respect to 
t. The constraint equation ([7]) now reads 

R- RS + f = 0, (14) 

where T = T vv = T uu is given by 

f = K[p-p 2 + z], (15) 

and z = z v = z u . Note that the equalities T vv = T uu 
and z v = z u are direct consequences of the homogeneous 
set-up. The latter equality also implies that z must be a 
constant. (By homogeneity z, could at most depend on t, 
but then dz/dt = z >v = z u>v = 0.) 
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The Schwarzschild solution inside the BH may be ex- 
pressed in the explicitly homogeneous form 



R 



M, 



S = t. 

In particular it satisfies the relations 
R= -e s = R- M, 



(16) 



(17) 



which will be useful in the analysis below. Note also the 
classical homogeneous constraint equation R — RS = 0. 



FIG. 2: A numerical plot of R (dashed) and of R (solid) near 
the R = K singularity, for K = 1 and M = 10 2 . While R 
is bounded, the divergence of R at the singularity is clearly 
seen. 



SINGULARITIES 

The semiclassical equations are singular at the points 
R = and R = K, as can be seen from Eqs. (| 121 1 3[) . In 
this section we will look at the asymptotic behavior near 
these two special R values. 

Since we set the initial conditions for the homogeneous 
solution at R » K, the first singularity to be encoun- 
tered is R = K . As it turns out, this singularity is char- 
acterized by the divergence of R and S, while R and S 
are finite. 

We therefore assume that in the right-hand side of Eqs. 
(|12I13|) we can neglect the terms proportional to e com- 
pared to those oc R 2 , obtaining at leading order 

R = - 2 (R-K) ' ( 18 ) 
& = 2K{R-K) ■ (19) 

Eq. (fT5|) constitutes a closed equation for R, which is 
easily solved: 

R(t) = K + B\t-t Q \i . (20) 
Then Eq. fT5]) is solved to yield 

S(t) = -§\t-to\~ 3 + (t-t )B 2 + B 3 . (21) 

The solution depends on four free parameters 
B, B 2 , B 3 ,to as required, hence this is a locally-generic 
asymptotic solution. We see that the variables R, S are 
continuous at R = K, but R diverges as \t — iol" 1 ^ 3 (and 
the same for S). This behavior, which is numerically 
demonstrated in Fig. [21 justifies our preassumption that 
e s is indeed negligible compare to R 2 . 

In the second singular point R = 0, we expand Eqs. 
(| 1 21 13[) near R = under the same assumption as before, 
namely that the terms proportional to e s may be ignored 
(this requires that e s << R 2 /R, which will again be 
justified a posteriori). We obtain 

i?.= i, (22) 
S = -£t. (23) 



These equations are easily solved: 

= d |t- i | 2 (24) 

and 

S (t) = C 2 + C 3 \t-t \ + 2 In \t - t , | (25) 

where Ci, C2, C3, to are four arbitrary parameters. We 
find that the present situation differs from the R = K 
singularity, because now R (like R) is finite and S (like S) 
diverges. Nevertheless, our preassumption e s « R 2 /R 
is still justified, because e 5 vanishes as (t — t ) 2 whereas 
R 2 /R approaches a constant. 

It is quite surprising to find that although the evolution 
equations are singular at R — 0, the evolving solution for 
both R and e s is regular there. In particular, the metric 
function p = (S — lni?)/2 is finite at R = 0. Note, 
however, that the dilaton (j> diverges there. 

Our homogeneous solution starts at R = R >> K, 
and R decreases until the R = K singularity is reached. 
The divergence of R and S at R = K results in an inabil- 
ity to provide a unique prediction for the evolution of the 
fields beyond that point. Therefore the R = singularity 
and its neighborhood are beyond our domain of predic- 
tion. In the rest of this paper we shall concentrate on 
the R = K singularity, as well as on the global behavior 
in the domain K < R < M. 



CONSTRUCTING THE INTERIOR 
HOMOGENEOUS SOLUTION 

So far, we dealt with the local structure of the fields 
near the singularity. Our ultimate goal, however, is to 
understand the global behavior of the fields inside the 
BH. To this end, we need to follow the evolution of R(t) 
and S(t) from the initial hypersurface (say at R M) 
up to the singularity at R = K. We shall do this by de- 
signing a couple of analytical approximations, the union 
of which cover the entire domain K < R < M, We 
shall also augment the analytic approximations by di- 
rect numerical integration of the homogeneous evolution 



6 



equations (|12ll3p . Wc shall start by specifying the ini- 
tial conditions, which are required for both the numerical 
and analytical solutions. 



Setting the initial conditions 

The required initial conditions for the homogeneous so- 
lution are the four functions R, R, S, S, all set at a certain 
initial moment to- The basic strategy of setting the ini- 
tial data follows from the discussion in section,, Suppose 
that we want to explore the interior of the evaporating 
BH at a stage (i.e. v value) where its remaining mass 
is M (with M » K). Then we take the initial condi- 
tions to be those corresponding to the classical internal 
Schwarzschild solution with the same mass parameter M, 
at a hypersurface R = const = Rq, at the limit Rq — > M. 
From Eq. ([16]) this amounts to setting R = M, R = 0, 
S = 1, and S — to — > — oo. In actual numerical im- 
plementations we pick a large negative value of to (say, 
to = —30) and set S = to- Note that in such a horizon- 
limit setup of initial conditions the choice R = is cru- 
cial for regularity. Any other choice of initial R = at 
t — > — oo will lead to a solution which lacks a regular 
horizon, and is hence inappropriate for approximating 
the spacetime of an evaporating BH. 14 1 

For this initial-value setup, the constraint equation 
(|14[) implies that T vanishes at the horizon limit t — > — oo. 
Eq. (US]) then yields z = 1/4. 



Approximate solutions in the different regimes 

We need to evolve the fields R, S from the initial hy- 
persurface, which we set to be at R — > M, up to the 
singularity at R = K (recall we assume M >> K). We 
notice three important domains in this overall range : (i) 
the macroscopic domain R >> K. (ii) very close to the 
R = K singularity, and (iii) the domain K < R « M 
[Note that (iii) is an extension of (ii), and it overlaps 
with (i) .] We shall now discuss the approximate solution 
in each of these domains. 



The macroscopic (or classical) domain 

It is easy to see that in the domain R > > K the quan- 
tum contribution to the right-hand side of Eq. (|12[) is 
negligible compared to the classical contribution (K = 
0). The same applies to Eq. ([T3]). Therefore, in this do- 
main the semiclassical solution may be well approximated 
by the classical solution Q16|) . (For numerical verification 
see Fig. @] below.) 



Near the singularity 

This domain is characterized by the inequality R — 
K « K . The approximate solution in this region was 
constructed in the previous section, Eqs. ()20|21|) . This 
solution is characterized by four parameters, of which the 
most significant one is B. 

The small-R approximation (R « M ) 

We seek an intermediate approximate solution, in the 
domain where R is << M, yet it is not quite close to 
K (say, R ~ 3A'). Formulating the desired approxima- 
tion, that will apply both near the singularity and in the 
portion R << M of the classical domain, requires un- 
derstanding the common basis of these two approxima- 
tions. Consider the expression given in the right-hand 
side of Eq. (fT2|) for R. In the near-singularity approx- 
imation the first term in the right-hand side (which is 
oc e s ) is negligible compared to the second one since R di- 
verges, whereas S(t) is bounded as seen in Eq. (|2"Tj) . One 
observes, however, that this term remains unimportant 
even in the portion R « M of the classical domain. In 
this regime, both R and —e s are given approximately by 
R -M = —M, sec Eq. ([TTJ). hence e s /R 2 scales as 1/M. 
Therefore, if we fix R and increase M, the first term in 
the right-hand side of Eq. (JTSJ) becomes negligible — just 
like in the near-singularity approximation. Omitting this 
term wc obtain: 



R = -R' 



K 



2R (R - K) ' 



(26) 



We shall refer to this approximation as the small-R ap- 
proximation. Note that the near-singularity approxima- 
tion may be obtained from the small- R approximation by 
substituting R = K, as may be seen by comparing Eqs. 
(EG]) and OH). 

Solving Eq. [26] for R we obtain 



R = A 



R 



R- K' 



(27) 



where A is a free parameter. One more simple integra- 
tion yields an expression for t(r), which we do not need 
however. 

In a particular evaporating-BH space-time, the param- 
eter A in Eq. ([27]) is to be determined from the effective 
mass parameter M. This can be done by matching the 
small- R approximation and the classical approximation 
in their overlap domain K « R « M. Eqs. ([27]) 
and ([T7]) respectively yield R = A and R = —M in this 
domain, implying 



.4 



-M. 



(28) 



The small- R approximation was established above by 
justifying it in both edges of its domain of validity 
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(namely at R = K and at K « R « M). We still 
need to demonstrate its validity in the region in between, 
where R is of order a few times K. To this end we note 
from Eqs. (|27I28|) . that for fixed R and K, R scales as 
M (just like e s ). Therefore, in the right-hand side of Eq. 
(TT2"]) the second term scales as M 2 , whereas the first term 
scales as M and can therefore be neglected at the large- M 
limit (M >> K). This justifies the small- R approxima- 
tion a posteriori, in its entire domain K < R << M. 

The validity of the small- R approximation, and its ef- 
fectiveness compared to the near-singularity approxima- 
tion, is demonstrated numerically in Fig. [3J 




34.50 34.51 34.52 34.53 34.54 34.55 34.56 
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FIG. 3: Plots of R{t) for the full numerical solution (solid 
black) , the near-singularity approximation (dotted black) and 
the small- 7? approximation (dashed gray). All plots are for 
K = 1 and M = 10 2 . (a) The domain where K < R < 8K. 
The large deviation of the near-singularity approximation is 
apparent, (b) A zoom at the domain closer to the singularity, 
where K < R< 1.4K. 

Finally, we can now determine the near-singularity 
coefficient B in Eq. (|2H)) . as a function of the effec- 
tive mass M. Differentiating Eq. (f20|) yields R = 
-fB^VVl-R- !<)■ Substituting Eq. (H and R ^ 
K « M in Eq. (gT]) yields R = -MyjK/{R- K). 
Matching the two last expressions for R, we obtain 




This analytic expression for B (valid at the leading or- 
der in K/M) may be verified numerically by evaluating 
R^/R — K at the limit where R — > K. For K/M values of 
0.01, 0.003, and 0.001, we numerically find the fractional 
deviation of numerical B from (|29[) to be 0.016, 0.006, 



and 0.002, respectively — suggesting that this fractional 
deviation may scale as K/M. 

Global approximate solution 

A global approximate homogenous solution in the 
range K < R < M can be derived by merging the classi- 
cal approximation R = R — M and the small- R approx- 
imation R = — ^J~-r-tk t° obtain a single approximate 
expression: 

R={R-M)^^^. (30) 

Substituting R >> K in Eq. (|3"D|) yields the classical ap- 
proximation. Substituting R « M yields the small- R 
approximation (|2"T|) (recall A = —M). Since the union 
of the domains R >> K and R << M covers the en- 
tire range of integration K < R < M (recall that we 
assume M >> K throughout), the approximation ((30|) 
is valid everywhere. We shall refer to it as the large-mass 
approximation. 

A comparison of the large-mass approximation and the 
full numerical solution, in the entire domain R < M , is 
presented in Fig. [U 
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FIG. 4: Plots of R (a) and R (b) for the full numerical solu- 
tion (dashed black), for the global large-mass approximation 
(solid gray) , and for the classical Schwarzschild solution (dot- 
ted black), for K = 1 and M = 10 3 . The breakdown of the 
classical approximation for R at small R is evident. The va- 
lidity of the large-mass approximation throughout the entire 
domain K < R < M is clearly demonstrated. 
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EFFECTIVE LAGRANGIAN AND 
HAMILTONIAN FOR THE 
NEAR-SINGULARITY REGION 

The occurrence of a spacclike singularity at R = K 
makes it difficult (if not impossible) to explore the evolu- 
tion of physics — and spacetime — beyond R = K by semi- 
classical methods. In Ref. [f| it was suggested that quan- 
tum evolution will still be regular even at the (would-be) 
semiclassical singularity. This motivates one to study the 
evolution of R and S in the neighborhood of the singu- 
larity in a quantum-mechanical framework. One of the 
simplest options would be to address the quantum prob- 
lem at the " mini-superspace" level: Since at the semi- 
classical level the homogeneous solution provides a good 
approximation to the actual spacetime evolution (as long 
as M » K), one may drop all spatial derivatives and 
analyze the simpler problem in which R and S depend on 
t solely. Furthermore, since the critical stage of evolution 
is the transition across the R = K singularity, it may be 
sufficient (at least as a first step) to analyze the problem 
at the leading order in R — K. 

This problem of quantum evolution is certainly beyond 
the scope of the present paper. Nevertheless, we do take 
here the preparatory steps which are to be implemented 
already at the semiclassical level: Namely, the formula- 
tion of an effective Lagrangian (and subsequently an ef- 
fective Hamiltonian) that describes the evolution of R(t) 
and S(t) near R — K. We therefore seek a Lagrangian 
L(R, S, R, S) that will recover the evolution equations 
(fTglTgj) near R = K. 

Since Eq. (JT8J) constitutes a closed equation of motion 
for R(t) , we start by constructing an effective Lagrangian 
L(R, R) for this equation. This turns out to be 

L R = ±(R-K)R 2 , (31) 

as one can easily verify by applying to it the Euler- 
Lagrangc equation. The momentum conjugate to R is 

P R = {R- K)R. (32) 

This yields the effective Hamiltonian for the i?-ficld: 



0- 

Next we treat the combined system of R and S (still 
in the homogeneous framework and at the leading order 
in R — K) . To this end we find it convenient to replace 
S by the new variable Z = R + KS. From Eqs. (|18I19|) 
it follows that Z = 0. Hence Z(t) corresponds to a free 
one-dimensional motion, with Lz = \Z 1 , conjugate mo- 
mentum Pz = Z, and Hamiltonian Hz = \P"z- The 
overall near-singularity Lagrangian and Hamiltonian are 
L = Lr + L z and H = H R + H z . 



DISCUSSION 

In this manuscript we primarily discussed two closely 
related issues. The first is the structure and main prop- 
erties of the spacelike singularity which forms inside a 
CGHS [H evaporating BH. The second is the structure 
of the fields (spacetime and dilaton) in the entire range 
of the BH interior, from the horizon up to the singularity. 

As a first step we established the homogenous approx- 
imation, a central tool which allowed us to analyze the 
above issues by means of ordinary differential equations. 
This approximation is valid as long as the BH mass is 
macroscopic. 

We find that the field equations admit two singular 
points: at R = and R = K . Using the homogenous 
approximation, we wrote down the asymptotic form of 
the field equations near each of these singularities [Eqs. 
(| 1 81 19[) and (|22I23[) ] and solved them exactly, obtaining 
the asymptotic solutions (|20I21[) and (|24I25|) at the two 
singularities. In the scenario of BH formation, R starts at 
large values and then decreases, so the R = K singularity 
will be reached first. It is therefore not clear whether the 
R = singularity will ever form. Hence, in the rest of 
our analysis we concentrate on the R = K singularity. 

Next we explored the evolution of the fields inside the 
BH from the horizon up to the singularity. By combin- 
ing several approximations, we composed a global ap- 
proximation for R as a function of R, given in Eq. (|30[) . 
This expression can be further integrated in a straightfor- 
ward manner to yield an expression for R(t) (though in a 
functionally- implicit form). It should also be possible to 
use this expression in order to construct an approximate 
global solution for S(t), but this requires some technical 
work which is beyond the scope of this manuscript [ItJ ■ 

The overall approximate expression (|30|) for R(R) al- 
lows us to determine the parameter B of Eq. (|20j) , which 
was left as a free parameter by the local analysis near 
R = K. This parameter is found to depend on the BH's 
effective mass M, through Eq. (|2"5|) . 

Our analysis made an extensive use of the local ap- 
proximate homogeneity which characterizes the R = K 
singularity. One should bare in mind, however, that the 
overall, large-scale, structure of the singularity is inho- 
mogeneous after all. The parameter B, which character- 
izes the local singularity structure, slowly drifts along the 
R = K singularity line of Fig[TJ This is explicitly seen 
by substituting M = M(v) in Eq. (|29j), (A similar drift 
may also apply to the other local parameters B2, -B3, to.) 
This qualitative picture holds as long as the remaining 
BH mass is macroscopic, namely M(v) >> K. 

From the asymptotic solution (|20I21[) , one immediately 
observes that R and S are continuous at R = K, yet their 
derivatives diverge there. The same applies to p. Thus, 
the metric tensor is continuous and non-singular, but its 
derivative diverges. The R = K singularity may be clas- 
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sificd as deformationally-weak [H , [HI : An extended body 
will only experience a finite tidal deformation as it ap- 
proaches the singularity. 

We briefly discuss here the possibility of extending 
semiclassical spacetime beyond the R = K singularity. 
From the discussion above it is obvious that a differen- 
tiable extended spacetime is ruled out, because differen- 
tiability already breaks down on the approach to R = K 
from the past. On the other hand, continuous exten- 
sions do exist. In fact, there is a two-parameter family of 
such continuous extensions (characterized by the values 
of B and B2 beyond R = K). Thus, in our view, the 
problem of extending the semiclassical solution beyond 
the singularity is primarily a problem of ambiguity: A 
priori it is not clear which of the infinite possible exten- 
sions will be chosen by Nature (if any). We expect that 
a fully-quantized treatment (like [B| for example) should 
provide the answer to the extension problem: First, it 
should tell us if a semiclassical phase at all develops be- 
yond the R = K singularity, and second, if such a semi- 
classical extension indeed develops, it should determine 
which semiclassical branch is selected. 

Finally we constructed an effective Lagrangian and 
Hamiltonian near the R — K singularity. Such an ef- 
fective Hamiltonian may be a useful input for future 
quantum-mechanical treatments of the behavior near the 
singularity. 

This research was supported by the Israel Science 
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stant rate K/4 (at the leading order); hence, M(v) = 
Mi n u - (K/4)(v - Vinit), where M init and Vi n u respec- 
tively denote the initial mass and the moment v of col- 
lapse. 

[12] The main effect of the inhomogeneity is to slowly drift 
the parameters of the homogeneous solution [e.g. the pa- 
rameters B and to in Eq. (|20p ] in the spatial direction 
along the singular line. The local temporal structure is 
not affected at the leading order. This was verified analyt- 
ically by expanding the inhomogeneous singular solution 
off the moment of singularity to (x) (at which R = K) up 
to seven orders in t — to(x) [3 

[13] Note that the above argument which establishes the 
homogeneous approximation does not directly apply to 
the case Rq = Mo, because the latter hypersurface in 
Schwarzschild is null rather than spacelike. Nevertheless, 
the limiting solution obtained as Ro — * Mo is well defined 
and is perfectly regular at the horizon. 

[14] The homogeneous evolution equations (I12I13|I admit a 
generic class of solutions in which 7? approaches a nega- 
tive constant as t — > —00. However in this class obviously 
R — > 00 at that limit, and the curvature diverges as well. 

[15] This is the unique z = z v — z u value in a precisely ho- 
mogeneous solution with a regular horizon. We point out, 
however, that the actual evaporating-BH solution is not 
precisely homogenous, and one finds that at the horizon 
limit z u — 1/4 and z v = 0. This difference between z u 
and z v may be regarded as a measure for the deviation of 
the actual evaporating-BH spacetime from the precisely 
homogeneous solution. Simple analytical arguments sug- 
gest that the effect of this change in z by 1/4 will be 
small for large M/K (and vanish at the macroscopic limit 
K/M — > 0). We explored this effect numerically (within 
the homogeneous framework), for K/M ranging between 
0.01 and 0.001, and verified that this is indeed the case. 
As a quantitative measure for this effect, one can look 
at the change in the actual B parameter at the R = K 
singularity, induced by such a change of 1/4 in z. Our 
numerical results suggest that the fractional change in B 
scales as K/M. 

[16] The effective Hamiltonian Ha may also be obtained in a 
more direct manner as follows: One starts from the orig- 
inal Lagrangian density fl]), omits all spatial derivatives 
(namely, replace d v and d u by a t-derivative), omit the 
fi contribution, and transform (j>, p to R, S, to obtain 
a Lagrangian L(R, S, R, S). Then one constructs from 
L the conjugate momenta Pr,Ps and the Hamiltonian 
H(R,S, Pr, Ps)- One then expands H in powers of R— K. 
The effective Hamiltonian Hr is obtained (up to a multi- 
plicative constant) as the leading order in this expansion. 

[17] The specific form of S(t) is less relevant here than that 
of R(t), since the singularity is determined by the value 
of R but insensitive to S. 



